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A consistent perturbation theory expansion is presented for 
phase-ordering kinetics in the case of a nonconserved scalar 
order parameter. At zeroth order in this expansion one ob- 
tains the theory due to Ohta, Jasnow and Kawasaki (OJK). 
At the next nontrivial order in the expansion, worked out in 
d dimensions, one has small corrections to the OJK result for 
the nonequilibrium exponent A and the introduction of a new 
exponent v governing the algebraic component of the decay of 
the order parameter scaling function at large scaled distances. 



I. INTRODUCTION 

Significant progress has been made on the theory of 
phase-ordering kineticsEJ using methods that introduce 
auxiliary fields that are taken to have gaussian statis- 
tics. These theories well describe the qualitative scaling 
features of ordering in unstable systems. The methods 
developed in OJKa and TUG0 each have appealing as- 
pects and separately give good descriptions of different 
aspects of the ordering problem. A major lingering ques- 
tion is why these methods work as well as they do and 
how they can he .reconciled and improved. Thus there has 
been a searchaH for the field theory description where 
OJK or TUG is the zeroth-order approximation in some 
systematic .expansion. Such an expansion is presented 
in this paper!! It has the OJK result as its zeroth-order 
approximation. More importantly it indicates how one 
goes forward to improve on these theories. The case of 
a scalar order parameter is treated in this paper. These 
ideas are generalized in a companion paper to the case 
of the n-vector model and systems with continuous sym- 
metry in the disordered state. 

The problem of interest is the restoration of equilib- 
rium in a system rendered unstable by a rapid temper- 
ature quench to a regime where the final state corre- 
sponds to a broken discrete symmetry. The ordering is 
controlled in this case by the decreasing area of domain 
walls separating the competing final degenerate states. 
The two coexisting theories, OJK and TUG, have been 
useful in understanding certain aspects of this problem. 
The theory developed in TUG has led to nontrivial ex- 
pressions for the nonequilibrium exponent A, defined in 
detail below, which are in good agreement with values 



known exactly or from simulation-data. However, as dis- 
cussed by Mazenko and WickhamQ, in the TUG approach 
the auxiliary-field correlation function exhibits a non- 
analytic structure at short-scaled distance which leads 
to unphysical jjeaults when used in calculations of— de-, 
feet densititesBliJ and defect velocity distributiona!3EHl 
for systems with continuous symmetry in the disordered 
states. More recently they showedt3 that the OJK the- 
ory can be derived from the exact continuity equation for 
the defect densities for point and line defects and leads 
to smooth physical results for defect properties. On the 
other hand, the OJK result is only compatible with rather 
trivial results for the exponent A. Thus, at this point, one 
does not have a theory which is both smooth enough for 
treating defect dynamics and yet robust enough to give 
nontrivial results for the nonequilibrium exponent A. 

In principle the task for the theorist in this problem 
may seem obvious: Linearize the order parameter equa- 
tion of motion and formally arrive at a gaussian field 
theory. Then do perturbation theory in the remaining 
nonlincarity. This is the well known path in conventional 
field theory. The problem, however, becomes clear when 
one looks at the simplest form for the equation of motion 
satisfied by the order parameter if> in dimensionless units 
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ip -ip 3 + V 2 ip 



(1) 



It is clear from this form that there is no dimensionless 
nonlinear coupling in which to expand. If one expands 
in the nonlinear term ip 3 one obtains, exponential growth 
in time (the Cahn-Hilliard theorylij in the case where 
one has a conserved order parameter) which, as a zeroth- 
order approximation, does not include any of the basic 
qualitative features of the long-time ordering. The rea- 
son is that the nonlinearity is essential in stabilizing the 
growth in the long-time limit. Indeed the combination 
ip — ip 3 must become small as the system orders. A mean 
field theory by, Langer, Bar-on, and Millenla, is an im- 
provement but is ultimately flawed by the inability to 
treat the separation of the two characteristic lengths in 
the problem. The dominant scaling length, L(t), grows 
algebraically with time and characterizes the average sep- 
aration between defects in the system. The other length 
is the equilibrium correlation length £g. Clearly at suf- 
ficiently long times L(t) ^> £g. This two length-problem 
was addressed by Mazenko, Vails and Zannettilla (MVZ) 
who argued that the solution to this problem is to sepa- 
rate the order parameter field ip into a peak contribution 
a and a fluctuating contribution u, ip = a + u . Then a 
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is associated with ordering on the length scale L{t) and 
u with equilibration on the length scale £g. In the struc- 
ture factor the a variable is identified with the growth 
of a Bragg peak with width and u is identified 

with the Ornstein-Zernike contribution with width gov- 
erned by Cij 1 - One reason for going over to the auxiliary 
field method is to insure that the Bragg peak grows with 
the proper weight which corresponds to the equilibrium 
average of the order parameter squared. 

The work of OJK predates that of MVZ but fits into 
the picture developed there. OJK ignore the fluctuations 
(it = 0) and assume a is a function of a gaussian auxil- 
iary field which is diffusive inJiatU|re. There have been a 
number of subsequent papera!3il3 clarifying the nature 
of the OJK result. More recently, in the work of Bray 
and HumayunB, there has been an effort to derive the 
OJK theory in a systematic fashion. This approach will 
be discussed below in section III. 

An alternative implementation of the ideas of MVZ 
was carried out in TUG. In this discussed below, 

the mapping of the order parameter onto an auxiliary 
field m is motivated by the idea that m measures the 
distance to the nearest defect. Coupled with the assump- 
tion that the auxiliary field is gaussian, the equation of 
motion satisfied by a is enforced, and one obtains a the- 
ory describing most of the features satisfied by the order 
parameter scaling function. 

All of these theories with auxiliary fields satisifying 
gaussian statistics have beenJiimped together as gaus- 
sian closure approximations^. It has been clear for 
some time -that we have needed to have theories which 
go furthcrEjE3|-||-, Previous efforts at a post-gaussian 
approximationstffl had some success but are difficult to 
control. 

How does one construct a perturbation theory expan- 
sion for the auxiliary field m? This theory will be unusual 
since, on average, m is growing with time. Thus standard 
polynomial nonlinearities in m would be a problem. As 
discussed in section III there are some difficult technical 
problems in developing on expansion method which is 
self-consistent. The resolution to this problem is to first 
choose the proper introduction of the auxiliary field m, 
and then organize the treatment of the associated non- 
standard nonlinear field theory through the introdution 
of an expansion parameter, <fi p . At zeroth order in this 
parameter one obtains the theory due to OJK and at 
second order one finds expressions for the corrections to 
the nontrivial exponents characterizing the ordering in 
these systems. This parameter, (f> p , is associated with 
the unusual nature of the nonlinearity in the problem. 
Instead of a polynomial nonlinearity one has the sign of 
the field giving the driving nonlinear terms in the equa- 
tion of motion for the auxiliary field. This expansion 
in 4> p appears well behaved order by order in perturba- 
tion theory and in some ways is similar to the bare ex- 
pansion in the quartic coupling in </> 4 field theory. One 
has, for example, the possibility of resummation or use of 
renormalization group methods. We appear to be fortu- 



nate in this case since lower-order approximations appear 
to work well, give reasonable results for anamolous di- 
mensions, without the need for extensive reorganization. 
One encounters a rather straightforward exponentiation 
of logarithmic divergences. Unlike in critical phenomena 
the logs appear for all dimensionality d. In this case they 
are driven by internal time integrations. 

A key assertion in the theories developed previously is 
that the auxiliary field can be treated as having gaussian 
statistics. It should be understood that a theory with this 
feature at zeroth order must have the property that all of 
the higher order cumulants for the field m must vanish at 
this order. It is shown here that indeed the n-point cu- 
mulant is of — I) in the expansion parameter. Thus 
the two-point cumulant, G2, which enters the determina- 
tion of the order-parameter correlation function at lowest 
order, is of O(0) as expected. The four-point cumulant 
is of 0(1) and so on. This has the consequence that 
any function of the auxiliary field m can be expressed in 
terms of these cumulants and evaluated in perturbation 
theory 

Turning to the question of the existence of a small 
parameter, it appears that the most direct connection 
of the expansion in <fi p to a more conventional expan- 
sion parameter is to the 1/n expansion in the n— vector 
model. This will be discussed in detail in paper II in 
this series. This expansion does not appear to be related 
to a 1/d expansions. Thus the method developed here 
has more in common with methods in critical phenom- 
ena where one works in the spatial dimension of interest 
as compared, for example, to the e-expansion about four 
dimensions. 



II. OVERVIEW 

A. Setting Up the Problem 

The system studied here is the domain-wall dynam- 
ics generated by the time-dependent Ginzburg-Landau 
(TDGL) model satisfied by a nonconserved scalar order 
parameter ?/>(r, i): 



(2) 



where T is a kinetic coefficient, F is a Ginzburg-Landau 
effective free energy assumed to be of the form 



d d r[ ^(VV0 2 + V(V>) 



(3) 



where c > and the potential V is assumed to be of the 
symmetric degenerate double-well form. We expect only 
these general properties of V will be important in what 
follows. r\ is a thermal noise which is related to T by 
a fluctuation-dissipation theorem. We assume that the 
quench is from a high temperature (Tj > T c ), where the 
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system is disordered, to zero temperature-where the noise 
can be set to zero (77 = 0). It is believectl that our final 
results are independent of the exact nature of the initial 
state, provided it is a disordered state with short-ranged 
correlations. 

If we rescale lengths and times we can put our equation 
of motion in the dimcnsionlcss form 



A(1)V(1) = -V'[V(1)] 
where the diffusion operator 



(4) 



(5) 



is introduced along with the short-hand notation that 1 

denotes (r 1 ,ti). Eq.(nl) is just the special case of the 
— \r _ i„/,2 n,i,4 



potential V 



B. Summary 



It is well established!!] that for late times following a 
quench from the disordered to the ordered phase the dy- 
namics obey scaling and the system can be described in 
terms of a single growing length, L(t), which is character- 
istic of the spacing between defects. In this scaling regime 
the order-parameter correlation function has a universal 
scaling form 



C(12) = W1M2)> = ^(x,h/t 2 



(0) 



where ipo is the magnitude of the order-parameter in the 
ordered phase. The scaled length x is defined as x — (ri — 
7*2) /L(T) where, for the non-conserved order parameter 
case considered here, the growth law goes as L(T) ~ T 1//2 
where T = |(ti + I n the <case of the autocorrelation 
function f\ — r 2 — r we haveEffij 



T 



(7) 



where A is a nontrivial noncquilibrium exponent and ei- 
ther t\ or t 2 is much larger than the other. At equal-times 
and short-scaled distances 



T{x) = F(x,l) = 1- a|ar| + 



(8) 



This nonanalytic behavior as a function of x is indicative 
of Porod's lawEZI as conventionally given in terms of the 
Fourier transform 



HQ) * Q 



(9) 



for large scaled wavenumber Q. That all of the higher or- 
der terms-in Eq.(||) are odd in \x\ is known as the Tomita 
sum ruleES. It appears, as we will discuss below, that the 
large x behavior can, with proper definition of x, be put 
in the form 



Hx) « — e~i x * 



(10) 



where v is a nontrivial subdominant exponential. 

The main results determined in this paper is an explicit 
determination of T(x,t\jt2) in perturbation theory. At 
zeroth order we obtain the OJK result 



T{x,t 1 /t 2 ) = -sin 1 $o(*i,*2)e * 



where 



*o(*i,* 2 ) = 



T 



An 



(11) 



(12) 



where Ao = d/2 and = 0. Going to next order, 0(1), 
in the expansion, discussed in detail in section X, we find 
no change in the indices A and v but quantitative changes 
in F(x,ti/t 2 ). At 0(2) both A and v are shifted. The 
exponent A is given at 0(2) by 



A 



d 
2 



■ to 



2 d M d 

3 d/2+l 



(13) 



where the dimensionality dependent quantities, w, Kd 
and Md, are determined by 



2w + u d 2 a K, 



3 d/2+] 



d 

= 1+ 2 



dz 



yd/2—1 



o [(l + z)(3-z)]^ 



and 



f\ z*/ 2 - 1 ir 2 (d/2) 

M d = / dz- -j = 

' [i + z} d 2 r(d) 



(14) 



(15) 



(16) 



The exponent v governing the algebraic component of 
the large x behavior is given at 0(2) by 



; 2 2 d+1 K d 



3 d/2+] 



(17) 



At lowest order one has the OJK results A = | and v 



0. 



While K d and oj can be worked out analytically for spe- 
cific values of d, the expressions are not very illuminating. 
Numerical values for A, v and uj are given in Tables I and 
II along with other results for comparison. 



III. REVIEW OF AUXILIARY FIELD METHODS 

How can one organize this problem in terms of a per- 
turbation theory expansion? Let us consider first the 
most direct method since the exercise is instructive and 
suggests other approaches. Let us ignore the fluctuation 
field u in the decomposition of the order-parameter field 
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and replace the order-parameter field ip with an auxiliary 
field m via the mapping 



(18) 



where cr[m], as in TUG, satisfies the Euler-Lagrange 
equation for the associated stationary interface problem 



d 2 a 
dm 2 



a 2 = V'[a[m}} 



(19) 



In this equation m is taken to be the coordinate. It will 
generally be useful to introduce the notation 



dm e 



(20) 



A key point in the introduction of the field m is that the 
zeros of m locate the zeros of the order parameter and 
give the positions of interfaces in the system. Locally the 
magnitude of m gives the distance to the nearest defect. 
As a system coarsens and the distance between defects 
increases, the typical value of m increases linearly with 

L(t)- n 

Inserting the mappingpj given by Eq.(|l8|) into the 
equation of motion for i\> given by Eq. (0) one finds, after 
using the chain-rule for differentiation, an equation for 
m(l): 

(7i(l)A(l)m(l) = -a 2 (l) [1 - (Vm(l)) 2 ] . (21) 

To get an idea of the nature of this equation take the 
special case of a ip A potential where one can solve for 
the mapping analytically and obtain the usual interfacial 
kink form 



o~[m] = tanh(m/V2) 



(22) 



It is easy to show for this particular potential, a 2 — 
— \f2tjo\ , and the equation of motion for m can be writ- 
ten as 



A(l)m(l) = \/2<t(to(1))(1 - (Vm(l)) 2 



(23) 



The left-hand side looks like the diffusion equation while 
the right-hand side has two types of nonlincaritics. The 
first nonlinearity is o~(m) which looks like sgn(m) in the 
bulk. The second nonlinearity is given by the (Vm) 2 , 
term. Working along these lines, Bray and Humayunu 
(BH) made the assumption that one can make a cleaver 
choice of the potential and replace — o 2 (m) j <j\(m) by m 
and then work with the equation of motion 



A(l)m(l) = m(l)(l - (Vm(l)) 2 ) 



(24) 



This polynomial form looks promising . If one makes the 
assumption 



(l-(Vm(l)) 2 ) 



L 2 (h 



(25) 



then the resulting linearized equation for m is given by 



A(l)m(l) 



L 2 (h 



-m(l) 



(26) 



This equation is formally the same at the one we will find 
later in the zeroth-order theory developed here and does 
corresponds to the OJK theory. 

It is not clear how to obtain systematic corrections 
to Eq.(26). Consider Eq.(|24]) from the point of view of 
dimensional analysis. Since m sa L, and A rj L~ 2 , one 
has that the left-hand side Eq.(||) is of C(L _1 ), while the 
right-hand side is of 0(L). The only way that the right- 
hand side of Eq.(p4]) can be of C(L _1 ) is if there is the 
additional constraint, given by Eq.(|2^), that is enforced 
at all orders. Formally this reduces to the statement that 
the nonlinear interaction 



V(l) = m(l)(l - (Vm) 2 ) - 



L 2 (h 



-m(l) 



(27) 



must lead to self-energy corrections which are of C(L _1 ). 
Due to internal pairings of the Vm terms in such an ex- 
pansion, this constraint can not be enforced order by 
order and one arrives back at the dimensional analy- 
sis argument discussed above. Thus the assumption 
given by Eq_(E5|) is not self-consistent without further 
developmentEiT 

One ends up concluding that the m(Vm) 2 ) nonlinear- 
ity is causing the technical problems and the OJK theory 
is not the zeroth-order solution of this problem as posed. 
Since there is reason to believe that the OJK theory is 
a good approximation to the original problem, then one 
might conclude that the role of the m(Vm) 2 ) nonlinearity 
is technical and not crucial and one should reorganize the 
calculation so that it plays a less prominent role. Indeed 
it is useful to turn the argument around and suggest that 
it is the o~[m] nonlinearity in Eq.(p3|) which is important 
in the equation of motion for m and its role should be 
emphasized. 

Let us back up a bit. Suppose, instead of the rather 
rigid mapping given by Eq. ( |l8| ) , we follow MVZ and write 



(28) 



where o~[m] is still the solution to the Euler-Lagrange 
equation Eq.(|l9|) and u[m] is to be determined. Let us 
substitute this mapping into the equation of motion for 
the order parameter, use the chain-rule as in leading to 
Eq. (p3[) , and obtain 



A(l)u(l) + <7i(l)A(l)m(l) 



-y'[ f 7(l)+ U (l)]+a 2 (l)(Vm(l)) 2 



(29) 



Notice that the perspective is different here. This can 
be regarded as an equation for the field u. We then have 
the freedom to assume that m is driven by an equation of 
the type given by Eq.(03h. However we now choose this 
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equation such that we can find a self-consistent expansion 
about the OJK result. To begin we assume this equation 
of motion is of the form 

A(l)m(l) = $(h)a(m(l)) . (30) 

Then this is just Eq.(p3|) with (l-(Vra) 2 ) replaced by the 
time dependent quantity £ (t) . From dimensional analysis 
m sa L, A sa L~ 2 , a w L° so £ « L~ l . These restrictions 
are very important. Using Eq.(|o|) in Eq.(p9|) for u leads 
to an equation of motion for the field u{\): 

A(l)«(l) - -V'[a(l) + u(l)] + (7 2 (l)(Vm(l)) 2 

-<7i(l)e(ti)<7(m(l)) . (31) 

The most important aspect of the solution of the last 
equation for u — u[m] is that 

limimi^ooulm] = . (32) 

Indeed as far as the universal properties are concerned 
this is almost all we need to proceed. To understand that 
we can construct a solution for u with this property let us 
consider the special case where we have a ip 4 potential. 
The equation for u is given then by 

Au + (3cr 2 - l)u + 3cm 2 + u 3 

= -<T2 [1 - (Vm) 2 ] - o-&o . (33) 

In the limit of large \m\ the derivatives of a go expo- 
nentially to zero and the right hand side of Eq.(|33|) is 
exponentially small. Clearly we can construct a solution 
for u where it is small and linearize the left hand side. 
Remembering that a 2 = 1 away from interfaces in the 
bulk we have 

Au + 2u = -<t 2 (1 - (Vm) 2 ) - oi£ sgn(m) . (34) 

Notice on the left-hand side that u has acquired a mass(— 
2) and in the long-time long-distance limit the term 
where u is multiplied by a constant dominates the deriva- 
tive terms: 

2u = -cr 2 (l - (Vm) 2 ) . (35) 

We have dropped the term proportional to £(t) since it 
vanishes more rapidly than the other terms at large times. 
That the u field picks up a mass in the scaling limit can 
easily be seen to be a general feature of a wide class of 
potentials where c/q = V"[a = ±ipo] > . We have then 
on rather general principles that the field u must vanish 
rapidly as one moves into the bulk away from interfaces. 

One expects that the explicit construction of u is rather 
involved and depends on the details of the potential cho- 
sen. If we restrict our analysis to investigating univer- 
sal properties associated with bulk ordering we will not 



need to know the statistics of u explicitly. If we are inter- 
ested in determining interfacial properties then we need 
to know u in some detail. Thus for example if we want 
to determine the correlation function 

C^(12) =< (^ 2 (1) - ^ 2 )(V 2 (2) - ^ 2 ) > (36) 

we will need to know the statistics of the field u. However 
if we are interested in quantitites like 

C(12...n) =<<0(1)V>(2)---V(n) > (37) 

where the points 12. ..n are not constrained to be close 
together, then we do not need to know u in detail. Why 
is this? Consider, for example, 

C(12) =< V(1)V(2) >=< HI) + u(l)] [<r(2) + «(2)] > . 

The point is that because the field u{l) is nonzero only 
near interfaces, the average < u(l)er(2) > is down by a 
factor of 1/L 2 relative to < a(l)a(2) >. Out in the bulk 
we can use dimensional analysis to make the following 
estimates: (Vm) 2 w 0(1) , cr 2 w C(L" 2 ) , ai£{t)a w 
0(L~ 3 ) , Au w 0(L- 4 ) , where £ w L" 1 . Using these 
results we find that the equation of motion for u can be 
put into the form 

V'[a + u] = cr 2 (Vm) 2 . (38) 

The usefulness of this equation is that it allows us to 
express the equation of motion for the order parameter, 
in the bulk, in terms of the field m, 

AtP = -cr 2 (Vm) 2 . (39) 

We will demonstrate the usefulness of this result in sec- 
tion XI. Our picture therefore has the m field driving the 
order parameter and the m field satisfying the nonlinear 
Eq.©. 

IV. FIELD THEORY FOR AUXILIARY FIELD 

Let us consider the field theory associated with the 
equation of motion for m(r,t). Our development will 
follow the standard Martin-Siggia-RoseEa method in its 
functional-integral form as developed by DeDominicis 
and Pelittl. The analysis begins with the basic equa- 
tion of motion for the field m given by Eq. (^o|) . In the 
MSR method the field theoretical development requires a 
doubling of operators to include the field M which is con- 
jugate to to. We also organize things so that the initial 
field TOo(r) is also treated as an independent field. Thus 
it is assumed that m is zero for t < to and one must add 
a term 5(ti — to)nio(ri) to the right-hand side of Eq. (p^p|) . 

Following standard procedures, averages of interest are 
given as functional integrals over the fields m, M and mo 
weighted by an action A: 

< m(l)m(2) • • • m(n)M(n + l)M (n + 2) ■ ■ ■ M (n + I) > (40) 
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VmVMVm m(l)m(2) ■ ■ ■ m(n) 



x M(n + l)M{n + 2) • • • M(n + |) e ^(m,M, roo ) j z 
where 

Z = J VmVMVm e A ^ M ^ . (41) 
The action takes the form 



A{m,M,m ) = -i / dlM(l) 



x [A(l)m(l) - ^(l)a(l) - 5{h ~ t )m (l)} 



\i d d n [ d d r 2 m (f 1 ) 5 - 1 (r 1 -r 2 )mo(f 2 ) (42) 



where we use the notation, J dl = J dtid d ri , and where 
we assume, as is appropriate in this case, that the initial 
field is gaussian and has a variance given by 



< m (fi)m (f 2 ) >= g(n - r 2 ) 



(43) 



We will not have to be very specific about the form of the 
initial correlation function g. It will be very convenient 
to generate our correlation functions as functional deriva- 
tives in terms of sources which couple to the conjugate 
fields. Thus we introduce 



S[h, H] = exp J dl [ft(l)m(l) + H{1)M{\)\ (44) 
and define 

Z[h,H] = J VmVMVm e A{m - M - mo) S[h, H] 

= J VmVMVm Q e AT{mM ^ ma) (45) 
where the total action is defined 

A T = A + I dl [h(l)m(l) + H{\)M{\)\ . (46) 



The fundamental equations of motion are given by the 
identities 



/ 



VmVMVmo 



I 



VmVMVm 



VmVMVmc 



8M{1) 
5 

5m(l) ' 
S 

Sm (l) 



^A T (m,M,m ) _ g 



,AT(m,M,mo) _ q 



^A T {m,M,m ) 



which reduce to 



< 



8M(T) 



A T (m,M,m ) > h = 



< 



5m(l) 



< 



Sm (l) 



A T (m,M,m ) > h = 



A T (m,M,m ) > h = 



where the subscript h indicates that the average includes 
the source fields h and H: 



< ■■■ >h- 



j VmVMVm e AT ^ m ' M ' mo \ 
Z(h,H) 



Taking the functional derivative with respect to M just 
generates the original equation of motion Eq.(pO[) with 
an initial condition and a source term 

x 

v4^(m, M, mo) 



SM(1) 



= -i [A(l)m(l) - ^(l)a(l) - S(h - to)mo(l)] + H(l) . 
The functional derivative with respect to m is given by 

A r n 

■^-^A T (m, M, m Q ) = i [A(1)M (1) + ^1)<tm(1)\ + h(l 

where we have introduced the quantities 

d 



u M {l)=cx l {m{l))M{l) 
Finally we need the derivative 
5 



(47) 



Sm (l) 



Ax(m, M, mo) 



= iM(fi,t )- J d d r 2 g 1 (r 1 - f 2 )m {f 2 ) . 

Inserting the results of taking these derivatives into the 
averages we obtain our fundamental equations: 

- i \k{l) < M(l) > h +^(1) < <tm(1) >h\ = Hi) (48) 



i [A(l) < m(l) > h < <j(1) > h -8{tx - t ) < m (l) > h ] 

= H(1) . (49) 

and 



G 



i < M(f u t ) > h = / d d r 2 g 1 {ri - r 2 ) < m a (r 2 ) > h 



The last equation allows one to solve for the average of 
the initial field in terms of the MSR field M, 



j-ggjClearly we can go on and generate equations for all of 
the cumulants by taking functional derivatives. Let us 
introduce the notation that G J 4 1) ^ 2l .... l ^ n (12...n) is the 
order cumulant for the set of fields {Ai,A2, A n } 



, / h 



< m (fi) > h = i / d a r 2 g{f x - f 2 ) < M(r 2 ,t ) >h 



where field A\ has argument (1), field A 2 has argument 
(2), et cetera. This notation is needed when we mix 
(51) cumulants with m and M. As an example 



Eq.(|49|) can then be written in the form 

i[A(l)<m(l) > h -f (1) < <r(l) >/J 

= - J d2 n (12) < M (2) > h +i?(l) 

where 

n (12) = 6(h - t )<S(ti - i 2 )a(fi - fa) . 



(52) 



(53) 



All correlation functions of interest can be generated as 
functional derivatives of < m(l) >h or < M(l) >h with 
respect to and H(l). 

In the limit in which the source fields vanish, each term 
in the two fundamental equations vanish. Therefore it 
is derivatives of these equations which are of interest. 
Taking the functional derivative of Eq. ( |48| ) with respect 
to h(2) gives the equation for the response function 



GAf m (12) 



< M(l) > h 



with 



A(l)GW(12)+f(l) 



5h{2) 

5 < a M (l) > h 



5h(2) 



(54) 



5(12) . (55) 



Taking the functional derivative of Eq. ( p2| ) with respect 
to H(2) gives 



A(l)G mM (12)-£(l) 



8 < g(l) > h 
SH(2) 



I 



d3 n (13)G M M(32)+o(12) 



(56) 



Taking the functional derivative of Eq. ( |52| ) with respect 
to h(2) gives 



A(l)G mro (12)-e(l) 



8 < ct(1) > h 



5h{2) 



= - jd'S n (13)G Mm (32) 



(57) 



Note that Eqs.(p5[) and (|56|) are redundant because of the 
relation 



G m M(12) = (?Mm (21) 



(58) 



G 



Mi 



,(1234) 



^ 3 < m(4) > h 
8H(l)Sh(2)8h(3) 



(59) 



As a short hand for cumulants involving only m fields we 
write 



G„(12---n) 



rn- 1 



8h(n)5h(n - 1) • • • 8h(2) 



< m(l) >h 



(60) 



The equations governing the n" 1 order cumulants are 
given by 



A(l)G Mn 



(I2...n) + Q n (12...n) 







(61) 



and 



i[A(l)G„(12...n)-Q„(12...n)] 



= - / dl n (ll)G Mm ...m(12...n) 
The Q's are defined by 
Q„(12...n) = £(1)- 



(62) 



5h{n)8h(n- l)---8h(2) 



t (12...n) = e(l) 



Sh(n)Sh(n- 1) ■ ■ ■ 8h(2) 



< <tm(1) >h (63) 



<tr(l)> h . (64) 



With this notation the equations determining the two- 
point functions can be written as 



A(l)G Mm (12) + Q 2 (12) =5(12) 



(65) 



*[A(1)G 2 (12)-Q 2 (12)] 



dl n (ll)G Mm (12) (66) 



We see that GMm...m and G mm ... m are coupled. 

The point now is to show that there is a consistent per- 
turbation expansion where the higher-order cumulants 
are also of higher order in some ordering parameter. To 
get started we need to express Qi(12) and Qi(12) in 
terms of the fundamental cumulants of the m-field. The 
first step in this direction is to show that these quantities 
can be expressed in terms of the probability distribution 



Ph(x, 1) =< 5(x - m(l)) > h = 



< 5(x -m(l))S[h,H] > 
Z[h,H] 



(67) 
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where the averages without the subscript h are weighted 
by A not At and the dependence on the source fields is 
explicit. The average over a(m) can then be written 



<a(m{l)) > h = / dxa(x)P h (x, 1) 



and 



Qi(l) - / dx £(l)a(x) P h (x,l) 



(68) 



(69) 



Evaluation of QUI) is only slightly more involved. We 
have, using Eq.(fl7|), that 



< o-m(I) >h=< M(l)cri(m(l)) > h 



< o- 1 (m(l))S'[/i,-H'] > 



(70) 



Z[h,H] 5H(1) 



Z[h,H] SH(1) 



< M(l) > ft + 



/ dx ai(x)P h {xA) 



dx ai(x)P h (x, 1) 



SH(l) 



Using this result in Eq.(63) with n = 1 



Qi(l) = / e(l)(7i(a;) 



< M(l) > h + 



SH(1) 



Ph(xA). 



Then any perturbation theory expansion for Ph{x, 1) will 
lead immediately to an expansion for Qi(l) and Qi(l). 
We can then obtain Q n and Q n by functional differenti- 
ation. 



V. PERTURBATION THEORY EXPANSION 

The perturbation theory expansion for Ph(x,l) is 
straightforward. Using the integral representation for the 
(5-function we have 



Ph(xA) 



dk 
2^ 



— ik 



where 



e~ lfc;r $(fc,/j,l) 



${k,h,l) =< e H(1) >h 



(71) 



(72) 



and H(l) = ifcm(l) . The average of the exponential is 
precisely of the form which can be rewritten in terms of 
cumulants: 



$(&;, h, 1) = exp 



OO -. 



(73) 



where G$ (1) is the s th order cumulant for the field 7Y(1). 
Since H(l) is proportional to m(l) these are, up to factors 
of ik to the s th power, just the cumulants for the m field: 



G%(l) = (ik) s G s (ll 
We can therefore write 

" (ik) 



$(k,h,l) 



exp 



,s=l 



G s (ll...l) 



(74) 



(75) 



Consider first the lowest-order contribution to Q n 
which does not vanish with the external fields h, H: 



Q 2 (12) = J dx £(l)a(x)-^-P h (x,l 



(76) 



We will assume, as we will show self-consistently, that in 
zero external field, n th order cumulants are of order § — 1 
in an expansion parameter we will develop. Expanding 
$(fc,/i, 1) in powers of the cumulants with n > 2 and 
keeping terms up to the 4-point cumulant, we obtain 



xPf(z,l) 



where 



and 



M) = / ^{k,hA)e-^ 



$ Q (k,h,l) = e ikG ^e-i k2G *( 11 



(77) 



(78) 



(79) 



Then, after taking the derivative with respect to h(2), 
setting the external fields to zero, and neglecting all cu- 
mulants with n > 2, we obtain 



$ (fc,/ l = 0,l) = e-5 fe2 ^( 1 ) , 



(80) 



and 



dk 



Q^ 0) (12) = y dxt{l)o{x) J ^-e- tkx ikG 2 (12)e-i k ' 2s ^ 

where we have defined in zero external field 

S a (l) = G 2 (U) =< m 2 (l) > . (81) 

There are several points to make here. Let us begin with 
the separation of the factor £(l)a(x) into a piece which 
contributes to the bulk universal properties and a part 
which does not. Note that in general we can write <r(x) 
in the form 



cr(x) — ipo sgn(x) + a(x) 



(82) 
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where ipo is the bulk ordered value of the magnitude of 
the order parameter and cr(x) goes to zero exponentially 
as la; | — > oo. For the ip 4 potential ipo = 1 and 



a(x) = —sgn(x) 



e V2\x\ + 1 



(83) 



This means that Q^\l2) can be written as the sum of 
two pieces: 



f)(12)=Q^(12) + Qr V; (12) 



)(0,B) 



where 



Q 2 °' B) (12)=£(l)V>oG 2 (12) 



J dx sgn(x) J ^-tke-^e-^ 3 ^ 



and 



(84) 

(85) 
(86) 

ifcz -Afc 2 S 2 (l) 



The key observation is that Q^' N \l2) is down by a fac- 
tor of l/L 2 relative to Q^' B \l2). It should be clear that 
this is a general result which will hold order by order in 
perturbation theory. As far as the bulk ordering proper- 
ties are concerned we can replace a(x) — > ipo sgn{x) . 
Thus we could have started with the equation of motion 
for the field m 

A(l)m(l) = £(l)Vo sgn(m(l)) (90) 

if we focus only on bulk ordering properties^. 

Turning to the other nonlinear quantity Q 2 , we have 
in general 



Q 2 (12) = / dx£(l)ai{x) 



G Mm (l2)+ < M(l) > h 



Sh(2) Sh(2)SH{l) 



(91) 



P h {x,l)- 



Q { °' N \l2) = G 2 (12) J dx £(l)flr(a;) / ^ike~ ikx e 
In both contributions we have the integral 



—e- ikx ike~* k3SaW = — f — e - ikx e~^ k2s2{1) 

2ir dx I 2ir 



Clearly term by term in our expansion we will find that 
the leading contribution comes from ui{x) — > ipo28{x) 
with the remaining terms leading to contributions which 
are of higher order in l/L. We therefore need only con- 
sider 

Qa(12)=£(lM) / dx25{x) 



dx 



*o(*M) 



(87) 



where $o(x, 1) is the Fourier transform of $o(k, h = 0, 1) 
dk 



G Mm (l2)+ < M(l) > h 



<5/i(2) Sh(2)8H(l) 



P h {x,l) 



$ (x,l)= / * e - lfca: e-^- 2 s 2 (i) 

J 27T 



e 2s 2 (i) 

V 2 ^2(l) 



In evaluating Q^' N \l2) we can take the derivative of 
&o(x, 1) and expand in inverse powers of 5^(1) to obtain 
to leading order 

Q^ N \l2) = G 2 (12) / dx^l)a{x) I 

J y2irS 2 (1) 

where it is crucial that cr(x) vanish for large |x| so that 
the x-integral exists. Turning to Q^ (12) we have 



if we are only interested in the bulk universal properties. 

A key observation is that as we analyze contributions 
to Q n or Q n we will find that each term consists of prod- 
ucts of correlation functions and response functions with 
legs tied together by factors defined by 



4> P {1) = I dx sgn(x) 



—ik 2 P +1 e 
2tt 



ikx 



*o(M) (92) 



Mi) 



Q 2 °' B) (12) =£(1)V G 2 (12) / dx sgn(x) 



dx 



$ (x,l) 



where we have used an integration by parts in going from 
the first to the second line and defined 



In this case we integrate by parts in the integral over x 
and use -£;sgn(x) = 25(x) to obtain 

g(°' B) (12)=e(l)V'oG 2 (12)2$o(0,l) 



<Mi) 



dk 

to* 



nS 2 (l) 



(93) 



£(1)^G 2 (12)4 



7r5 2 (l) 



(89) 



Each term in the perturbation theory expansion for Q n 
or Q n will be proportional to factors of <p p . The perturba- 
tion expansion is ordered by the sum of the labels p on <p p . 
Thus a contribution with insertions (fiifafii 1 each factor 
typically associated with different times, is of 0(4). We 
shall refer to this expansion as the (p- expansion. It should 
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be emphasized that at this stage that this is a formal ex- 



It should be clear that the vertices V s (l) are of at least 



pansion. At order n it is true that <p p ~ L ( 2 p +1 ) which O(s) in the (^-expansion. By direct expansion of ft 
is small, however it will be multiplied, depending on the 0, 1) about &o(k, h = 0, 1) we obtain 



quantity expanded, by positive factors of L(t) such that 
each term in the expansion in <f) p has the same overall 
leading power with respect to L(t). 

To see how this expansion works let us consider first 
the two-point quantity Q2(12), defined by 



v s (i)^(i) + ^4 +2 (i) Se{1) 



G! 



>s+3( 



Qa(12) 



dx £(l)sgn(x) 



SP h (x,l) 
Sh(2) 



rs2(l) , 5 8 (1) 



(94) 



2(4!) 2 



(101) 



Then since we will find SV(1) ~ C(| — 1), the terms in 
Using Eqs.@ and @ and taking the derivatives with the expansion for V s are of 0(s), 0(s + 3), 0(s + 5) and 



respect to ft. (2), we find, in the case of zero external fields, 

dk 



,(12) =^(1) / dx sgn(x) 



-ikx 



2tt 



ES^TTG 2s+2 (11...12) 



s=0 



(2s + l)\ 



${k,h = 0,1) 



(95) 



0(s + 6) respectively. 

Let us turn next to Q2(12) given by Eq. (|9l|) . In the 
limit of zero external fields we can set the term propor- 
tional to < M(l) > to zero and then evaluate the second 
derivatives, 5 2 Ph(x, l)/ffi(2\iff(l). After a significant 
amount of algebra we obtainL3: 



Qa(12) 



dk 

dx2ai)ip 6(x) I — e -^$(fc, 1) 

Z7T 



Since all odd cumulants vanish in the case of zero external 
fields and 



$(fc,/i = 0,l) 
where 



exp 



E 



(~i) s fc 2 



-5 2s (l) 



S 2 .(l) = G 2 a(H...l) • 

Let us define the set of vertices 
dk 



(96) 



(97) 



G ! Mm(12) + E — j - G 



mm...mMm\ I 



(» fc )V(s+2) 



< 

2tt 



-ikx 



H>W — I dx sgn(x) 



2tt 



1 ^—ikx 



ft = 0, 1) , (98) 



which reduces, after following the same set of steps in 
reducing the original expression for <^ p , to 



E 

s=0 



Q 

| mm...mMm 



(11. ..112) 



/dk 
— k 2 ^(k,h = 0,1) 



s=0 



(2s)\ 



(99) 



Before going on to discuss the structure of the per- 
turbation theory at higher order let us made sur e the 
which is independent of position. Then the quantity theory is sensible at zeroth order where, from Eqs.(lOC) 



Q2(12), which appears in the equation of motion for and ( J102| ) 
G*2(12), is given in the form 



Q 2 (12) = £(l)^oE 



s=0 



(~1) S 

(2s + l)! 



Q 2 0) (12) = C(l)^o(l)G 2 (12) 



G 2s+ 2(H...12) 



Q 2 0) (12) = £(l>/W>o(l)GA/ m (12) 



(103) 
(104) 



x j dx sgn(x) 

(-iy 



dk 
— i 

2tt 



- <fex ifc 2s+1 $(fc,l) 



VI. ZEROTH ORDER THEORY FOR TWO 
POINT CORRELATION FUNCTIONS 

The equations of motion at zeroth order for the two- 
point correlation function is given by Eqs.( |6l|) a nd ( [62] ) 
where we have used the definition of V s (l) given by with Q2 and Q2 replaced by Eqs.(104) and dl03[). Thus 



£(DE 



^ (2s + 1)1 



V s (l)G 2s+2 (11...12) 



(100) 



Eq. (B8h in the last step. 



we have 
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A(l) +Wo (l) G$ m (12) = 5(12) 



(105) 



Multiplying from the left by G^^ and using the defini- 
tion of the inverse gives the symmetric form: 



*[A(l)-^ (l)]G( 0) (12) = -y dl n (ll)G^ ro (12) (106) 
where we have defined 

wo(i)=£(i)iMo(i) • ( 107 ) 

The first step in the solution to these equations is to 
Fourier transform over space. Taking the equation for 
the response function first, we obtain 



-.(<>) 



G<°> (12) 



dl I d2 G^ M (ll)G^ M (22)n (12). (117) 



Taking the Fourier transform and inserting the results 
for the propagators and IIo we obtain 

-,(0) 



G 2 °\q,ht 2 ) = 0(h - t )9(t 2 - to)e-« C^" 2 '") 



x R(t 1 ,t )R(t 2 ,t )g(q) 



(118) 



G i ° ) m (^ih) = 6(t 1 -t 2 ) . (108) : 



where g(q) is the Fourier transform of the initial corre- 
lation function. Henceforth we will suppress writing the 



step functions in the correlation function. Inverting the 
Fourier transform we obtain 



This first-order differential equation has the solution 

rti 



G Mm(l, hh) = -i0(t 2 - h)exp 



dr (q 2 - u> (t)) 



G { 2 °\r,t l t 2 ) = R{t l ,t )R{t 2 ,t ) 
(109) 



d d q 
(2^ 



- 2 q 2 T 



(119) 



where it is convenient to introduce 



= -i9{t 2 -t 1 )R{t 2l t l )e- q2{t2 - tl) 
and we have defined 

Taking the inverse Fourier transform, using 



dd 1 e *? e -*<*-*) = e ~ i(H - t2> 



T 



to 



(120) 



(110) 



While we are primarily interested in the long-time scal- 
ing properties of our system, we can retain some control 
over the influence of initial conditions and still be able 
to carry out the analysis analytically if wc introduce the 
initial condition 



9(q) = 9oe-^ )2 



(2ir) d 
we obtain 



[47T(tl - t 2 )} d /* 



(111) or 



9(r) = .90 



e -\ixiif 
{2nP) d / 2 



(121) 



(122) 



GaLMi^) = ~i0(t 2 - h)R{t 2 ,h) 



e 4(*2-*i) 

[47T(t 2 - h)]^ 



(1 12^ nser ^ n ^ ^ orm Eq.( |ll9| ) and doing the wavenum- 
ber integration we obtain 



It is straightforward to show explicitly the result we ex- 
pect from symmetry considerations: 



G [ 2 °\?Mt 2 ) 



(123) 



,(o) 

T mM 



(r,t 1 t 2 ) = -i9(t 1 -t 2 )R(ti,t 2 ) 



e 4(*i-*2) 
[Mh - t 2 )] d / 2 



(n3) - R{tlMR{t2M [2n(p g +\T)}^ 
In the long-time limit this reduces to 



e 2' 



Let us turn our attention to the correlation function. It 
is useful to introduce the inverse propagators 



G { 2 0) (r,ht 2 ) = R(t u to)R(t 2 ,t Q )g 



<4m(12) = -i A(l) + wo(l) 5(12) 



e -r 2 /ST 

{8nT) d / 2 



(124) 



(%£(l2)=i[A(l)-w (l)]5(12) 



(114) 



(115) 



Let us turn now to the quantity R(t\,t 2 ) defined by 
Eq.( |110D . We have assumed that w 1/L(t) and we 
have 0o ~ 1 so we can write for long times 



which allows one to write the equation for the correlation 
function in the form 



t c + t 



(125) 



dlG$£(lW(l2) = - / dl n (ll)G- ro (12). (116) ^^-^Z^^IZ 



.(0), 



-,(()) 



where uj is a constant we will determine and t c is a short- 
time cut-off which depends on details of the earlier-time 



gral 
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dTW (r) 



■J 1 2 

we obtain 



dr- 



R(h,t 2 ) 



t c + T 



t 2 +t c 



lo In 



h+t c 
t 2 + t c 



, (126) 



(127) 



where 



and 



,9o 



(135) 



Inserting this result back into Eq.(124) leads to the ex- 
pression for the correlation function 



a = 4/f . (136) 
This equation is solved in the Appendix with the result 



G 2 X>) (r,t 1 t 2 )=g(0) 



t + t c 



t + t c 



e -r 2 /8T 
{%TlT) d / 2 



(128) 



2d! 



a(l + d/2) 



[l + a(t-t )]- d/2 



If we are to have a self-consistent scaling equation then 
the autocorrelation function (r = 0), at large equal times 
t\ = t 2 = i, given by 



Si°\t) = g 



-d/2_ 



t 



to+t. 



1 



(87ri) d / 2 



.9o 



(to + t c ) 2 " (Stt)- 1 / 2 



(129) 



must have the form s!f\t) — Agt. Comparing we see the 
exponent lo must be given by 



u> 



2 V 2' 



and the amplitude by 



1 



.90 



(t +t c )^ (8tt)<V 2 



(130) 



(131) 



One question which arises is whether the time depen- 
dence in u>o(t) should be viewed as externally driven or 
sclf-consistently developed. Rephrasing this question: 
Does the time evolution for m(t) evolve out of the early- 
time instability? It is instructive to see that this follows 
in a natural fashion. Suppose instead of Eq.(125) we as- 
sume 



u (t) = 

which is consistent with 
w [t) 



s 2 (t) 



)(*) 



2S 2 (t) 



(132) 



(133) 



since <pp( t) = \/2/nS 2 ( t). I nserting this result for uJa(t) 
into Eq.( |llO| ) and Eq.( |ll9[ ) with r = 0, one finds at 
equal times ti = t 2 = t, a self-consistent equation for 
the zeroth-order quantity (t) : 



Si°\t) = exp 2Cj 



(h 



t s£\t) I [l + a(t-t 



id/2 



(134) 



2Q 



a(l + d/2) 



[l + a(t-t )] 



(137) 



and ui is still undetermined. There are several interesting 
results which follow. If we look at the long-time limit of 
u>o(t) we obtain 



-(l + a(t-t ) 



(138) 



a(l + d/2) 

If we compare with Eq. ( |125| ) , we regain Eq. ( |130| ) and find 

t c =--t . (139) 
a 

Thus ou r pro cedure is c onsi stent and the system governed 
by Eqs.(118) and Eqs.(132) is driven to grow as desired. 
Since S$ (t ) is now a known function of time we can 
use Eq.(134) to evaluate R(t, to) and obtain the complete 
solution to this zeroth-order problem over the entire time 
regime 



G { 2 °\rMt 2 ) 



,9o 



[27r0? 2 + 4T)] d / 



. e -ir 2 /(£ 2 +4T) 



)sf 2 °\t 2)i 



1/2 



S 2 °\to) S^'ito 



(0), 



(l+a(ti -t )) d/2 (l + a(t 2 -t )) 



d/2 



In the long -tim e limit this can be written, using the def- 
inition Eq.( |l39| ) for t c , in the form 



G< 0) (r, ht 2 ) = A Q [(t x + t c ) (t 2 + i c )]* (1+ ^ - 



-r 2 /8T 



T 2 

with 



T d/ 2 



(140) 



Ao = 



2lD 



l + d/2 



(141) 



Notice that t his result for A differs from the result 
given by Eq.(131) and does not depend on the initial 
conditions. This is evidence that the coefficient Ao is 
nonuniversal. Note, if we choose to enforce the condition 
< (Vm) 2 >= 1 at late times (see Eq.(|25])) we can fix the 
parameter u> at this order with the result 
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2 M ^ 



(142) 



G Mm (l2)+ < M(l) > h 



8h{2) Sh(2)SH(l) 



Ph{x,l) 



The general expression for the correlation function can 
be rewritten in the convenient form 



Gf (r,ht 2 ) = \/sf (t 1 )4° 3 (* 2 )$o(ti*2)e-^ 2 /^ +4T) (143) 
where, using Eq. ( |139| ) for I /a, 

d/2 



$o(*i*a) 



y/(ti+t e )(t2+t c ) 
T + t c + to 



(144) 



The nonequilibrium exponent is defined in the long-time 
limit by 



^ 0) (*i)^ 0) (* 2 ) 



T + t + t c 



and we obtain the OJK result 

d 
2 



A, 



(145) 



(146) 



We put the subscript to on A to indicate the exponent 
associated with G 2 (12). In general G 2 (12) and the order 
parameter correlation function C(12) need not share the 
same exponents. 

Looking at equal times we have that 



(147) 



where the scaled length is defined by x = r/At . fo(x) is 
just the well known OJK result for the scaled auxiliary 
correlation function. 



dx2£(l)5(x) 



+ GMm (12) 



Gm in hi in 

(1234)^(^,1) 



Ph(x,l) 



Gm™( 13) 



G Mm (14) 



5 4 



5h(4)Sh(3) 
S 2 

5h(2)5h(4) 
S 2 

5h{2)Sh(3) 



Ph(x,l) 



Ph(x,l) 



Sh{A)Sh{3)Sh{2)SH{l 



-P h (x,l) 



(148) 



to 0(1). Assuming, as we will show self-consistently, that 
G n « 0(5 — 1) and V„ w 0(n), one easily finds that the 
0(1) contributions to Q 4 and Q\ are given by 



and 



Q 4 1) (1234) = w„(l)G 4 (1234) 
wi(l)G 2 (14)G 2 (13)G 2 (12) 



?i 1) (1234) = Wo (l)G Mmmm (1234) 



(149) 
(150) 



VII. PERTURBATION THEORY FOR 
FOUR-POINT CUMULANTS 

In order to see how things go at higher order in per- 
turbation theory it is useful to look at the lowest nonzero 
approximation for the four-point cumulants which enters 
into the 0(2) correction for the two-point cumulant. In 
order to compute G4(1234) and GMmmm(1234) to 0(1) 
in perturbation theory, we see from Eqs.(|^) and ( p2] ) 
that we must evaluate 



Q 4 (1234) 



and 



S 3 



8h(4)Sh(3)Sh(2) 
Q 4 (1234) = 

s 2 



£(l)ipo dx sgn(x)P h {x, 1), 



rOa(12) 



Sh(A)Sh(3) 



5h{A)5h(3) 



-wi(l) 



G M m(12)G 2 (13)G 2 (14) + G Mm (13)G 2 (12)G 2 (14) 



G M m(14)G 2 (12)G 2 (13) + G 2 (14)G 



Mmmm 



(1123) 



G 2 (13)G 



Mmmm 



(1124) + G 2 (12)G 



Mmmm 



(1134) 



where we have introduced the notation 
= £(1)^(1) . 



(151) 



For future reference it is not difficult to work out the 
0(2) contributions to Q 4 and Q 4 given by 



?i 2) (1234) 



wi(l) 



G 4 (1134)G 2 (12) + G 4 (1124)G 2 (13) 



13 



G 4 (1123)G 2 (14) 



(152) 



G 4 (1234) — — G^ M (\l)Iia{\2)G Mmmm 



(2234) 



and 



QT (1234) 



+ G^ M (11)(-^ 1 (1))G 2 (12)G 2 (I3)G 2 (14) 



(157) 



Inserting the result for GMmmm given by Eq.( |155| ) and 
G 4 (1134)G M m(12) + G 4 (1124)G M m(l 1 3iing Eq.tfLVfy gives the result 



+ G 4 (1123)G M ™(14) 



G 4 (1234) = G^ 0) (ll)[-iwi(l)] 



G Mm (i2)G 2 (13)G 2 (14) 



•«i(l) 



+ G Mm (13)G 2 (12)G 2 (14) + G Mm (14)G 2 (12)G 2 (13) 



G 



M 



mmmm 



(1134)G 2 (12) + G 



M 



mmmm 



(1124)G 2 (13) 



G 



Mrnmmm 



(1123)G 2 (14) 



(153) 



Using Eqs.( |l5C| ) and (fill), we easily find the determin- 
ing equation for GMmmm at lowest nontrivial order: 



Mi, 



,(1234) 



+ M(1) 



G Mm (12)G 2 (13)G 2 (14) 



(154) 



+ G Mm (13)G 2 (12)G 2 (14) + G Afm (14)G 2 (12)G 2 (13) = 
This equation is easily integrated to give 

Gm in in in 

(1234)= / dlG^ ro (ll)[-i Wl (l)] (155) 



G Mm (T2)G 2 (T3)G 2 (14) + G Mm (13)G 2 (12)G 2 (14) 



+ G Mm (14)G 2 (l2)G 2 (l3) 



Inserting Eq.(149) into Eq.(p2J), we have 



G mM (11)G 4 (1234) + iwi(l)G a (12)G 2 (13)G 2 (14) 



= -n (ii)G 



Mmmm 



(1234) 



(156) 



where an integration over repeated barred indices here 
and below is assumed. This can be integrated up to give 



G^ 3 m (ll)Ha; 1 (l))G 2 (12)G 2 (l3)G 2 (14) 



(158) 



If we evaluate all of the two-point correlations in G 4 and 
GMmmm at lowest order we obtain our lowest-order ap- 
proximation for the 4-point quantities. Notice that G 4 is 
properly symmetric under interchange of any two of its 
labels. 



VIII. STRUCTURE OF PERTURBATION 
THEORY AT HIGHER ORDER 

In order to understand the structure of the pertur- 
bation theory one can rather easily see that we should 
consider two classes of contributions to 



,(12.. .n) 



Sh(n)5h(n- 1)---Sh(2) 
dk 



2tt 



'*(*:, M) 



(159) 



The first class corresponds to all 2n — 1 derivatives with 
respect to h acting on Q(k, h, 1) to give, in the zero ex- 
ternal field limit, a product of two-point correlation func- 
tions of the form 

Q 2n (12---,2n) 



= e(l)^oV n _i(l)G 2 (12)G 2 (13) • ■ ■ G 2 (l, 2n) . 

There is another set of contributions where all of the 
derivatives except the first acts on the factor multiplying 
<fr{k, h, 1) and gives a contribution 

Q 2n (12 • • • , 2n) = £(l)VoV (l)G 2 „(123 • • • , 2n). 

Since each of these terms is of the same order in (f> pi we 
easily find the proposed result G 2p w 4> p -\ ■ A very 
similar analysis follows for the case of Q 2n . 
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IX. TWO-POINT CORRELATION FUNCTION 
AT HIGHER ORDER 



Using the results for GMmmm and G4 at 0(1) leads to 
the 0(2) result for the correlation function 



Following the same procedures it is easy to find that 
the next order contribution to Q 2 and Q 2 is of 0(2) and 
given by 



>' J| (I2. = -^G4(1112) 

'1 '(12) = ^— -C?Mmmm(lll2) 



(160) 



(161) 



G 2 (12) = G 2 0) (12) + G 2 Z ' 1; (12 



(2,1), 



+ G 2 24) (21) + G 2 2 ' 2) (12) 



where 



Gf 1} (12) = Gf(ll)S^ m (12)G^ m (22) (167) 



,(2) 



(0) 



and 



Inserting Eqs.(104) and (161) into Eq. (pq) gives the equa- 
tion for the response function to second order: 

GM m (ll)GM™(12) +i^-GMmmm(ni2) = 6(12). 

Similarly, the equation determining the correlation func- 
tion to second order is given by 



G-^°(ll)G 2 (12)+i^ 



G 4 (1112) = -n (lT)G Mm (12) 



These equations can be integrated to give 
Gm™(12) = G^ m (12) 



+ 4(11) — 1^ GMmmm (1112) 



and 



(162) 



Gf 2) (12) = -Ga(ll)n (2) (12)Gi° ) m (22) 



(168) 



The self-energy is the same as for the response function 
and 



1 



n (2) (12) = ~[-H(l)]Gr J (12) 3 H^(2)] 



(169) 



We need to evaluate G 2 2,1 ^ and G^" 2 \ The integra- 
tions over space, in d-dimensions, are straightforward 
since they involve products of displaced gaussians. Af- 
ter rescaling the internal time integrations t\ — Ty\, 
h = Ty 2 , T = i(ti + t 2 ), we obtain 

G { 2 2 ' 1} (12) = ^S^T)S^j2 d - 1 u J H (t 1 t 2 )J 1 (x,t 1 /T,T) 
rti/T 



J x {xM/T,T) 



to/T 



dyi I dy 2 Ri(y 1 ,y 2 )e 

to/T 



G 2 (12) = -G^(ll)n (l2)G Mm (22) 



+ ^(11)^G 4 (1112) 



Inserting Eq.(162) for Gwm into Eq.(163) gives 



(163) 



#1(2/1,2/2) 



d/2-l d/2-1 
Vl V 2 

l(yi + y2)(3 yi ~y 2 -( yi -y 2 y)} d /2 



G 2 2 < 2) (12) = y/S {l)S (2)—u J 2 9 (t 1 t 2 )J 2 (x,t 1 ,t 2 ) 



G 2 (12) = G 2 0j (12) + G 2 u; (ll 



(0), 



-iwi(l) 



G 



Mr, 



+ ^L(11)^|®G 4( 1112) 



,(1112) 



(164) 



J 2 (x,t 1} t 2 ) 



ti/T 



dyi 



to/T 



t 2 /T 



to/T 



dy 2 R 2 (yi,y 2 )e * 



Using our 0(1) result for GMmmm gives the 0(2) result 
for the response function 

G Mm (12) = Gil(12) + Gil.(ll)41(12)Gil(22) , (165) 



#2(2/1,2/2) 



d/2-1 d/2-1 
Vl 2/ 2 



[(2/1 + 2/ 2 ) 2 (3- 2/1-2/2)]^ 

where we have chosen x 2 = r 2 /4T and 

3yi - 2/2 



where the lowest-order self-energy contribution is given 
by 



31(2/1,2/2) 



3yi - 2/2 - (2/1 - 2/2) 2 



and 



S£l(12) = |[-ia, 1 (l)]GW n (]2iG^'.12r|../.-,(2i; . 1 IWii 



52(2/1,2/2) 



3 - ?/i - 2/2 
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The first thing we should do with this result is look 
at the contribution at this order to the onsite equal-time 
t\=t 2 = t correlation function given by 



S 2 (t) = Si 0) (t) + 2G^(11) + G^(ll) 



(2,1), 



•y(2,2), 



where we have the integrals 



1 + 2'V t /i(0,M) + — Lu 2 J 2 (0,t,t) 



/•l rVi 

Ji(0,l,i)= / dyi / dy 2 Ri(yi,y 2 ) 

J t jt J t jt 



J 2 (0,t,t) = dyi dy 2 R 2 {yi,y 2 ) . 

Jta/t Jto/t 

The key point here is that Ji(0, l,t) and J 2 (0, t, t) are 

logarithmically divergent at t — > oo. One can show to 
logarithmic order that 



Ji(0,M) = K d ln(t/to) + - 



where 



and 



where 



l! dz ¥ 



y d/2-l 



Mo,t,t) 



3 d/2 



+ z){'i - z)] d / 2 



M d ln{t/t ) 



M d 



I 



1 z d ' 2 - 1 _ 1 T 2 {d/2) 
Z [l + z] d ~ 2 T(d) 



We have then that 



s 2 (i) = si°\t) 



1 + w 2 2 d K d + 



M d 



3 d/2+l 

and a simple exponentiation of this result gives 



ln(t/t ) + ■ ■ 



(! + •••) • 



For self-consistency we must determine u> at this order. 
Remembering that at lowest order 



si°\t) = A t 2 "- d / 2 , 



we require that 



S 2 (l) = A l-) {l + -..) = At 

which determines u) at this order 



2lu + uj 2 2 d [ K d 



M d 



1 



3 d/2+l j " ' 2 

Then, for example, for d = 2 

K 2 = \lnZ , M 2 = 1 



(170) 



and 



v/1 + 2 n 3 + 4/9 - 1 

w = v , - — — -f = 0.687687370 . . . 

In 3 + 2/9 

For large d analytical progress leads to the results 

2 1 

K ^dY d+ - 



M d 



/2tt 1 

~d¥ 



+ ■ 



and 



— 1(^-0 ■ 



A numerical determination of u) shows that it is approx- 
imately linear with d over the whole range of d. 

The t/to — > oo singularity in G 2 (12) at 0(2) can be 
regulated by turning our attention from G 2 (12) to the 
quantities ^{t\t 2 ) and f{x,t\/t 2 ) defined by 



G 2 (12) = y/S{\)S{2)*(tit2)f{x, h/t 2 ) (171) 
and the constraint f(x = 0,ti/t 2 ) = 1 . If we write 



and 



then 



G 2 (12) = G|, 0) (12) + AG 2 (12) 

5 2 (l) = 5f (1) + A5 2 (1) , 
G 2 (12) 



Ht 1 t 2 )f(x,t 1 /t2) 



&°Ht u t 2 )fo(x) 1 



V5 2 (l)5 2 (2) 



AS 2 (1) , AS 2 (2) 



(172) 



4 0) (D 



sf(2) 



+ 



AG 2 (12) 



We can then separate the contribution to the on-site cor- 
relation function ^(ti, £2) from the general x-dependence 
and write 

*(ti,t 2 ) =$o(ti,t 2 )(l+^ 2 2 d A$(t 1 ,t 2 )) (173) 
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f(x, h/h) = f (x) (1 + uj 2 2 d W{x, h/h)) (174) 
where 



regular quantity for t\ 3> t 2 or t 2 ^> t\. Turning to $2, 
however, one finds that it is logarithmically divergent in 
these limits: 



(175) 



A 2Md 1 

$ 2 = - rT77r In 



2 /•*i/ T /-i/i 
= o / ^ / dy 2 Ri(yi,y2) 

L Jto/T J to IT 



3 d/2 — y T 
This means to second order we have 



*(tl,*2) 



d/2 



^ rt 2 /T nyi 

2 y d 2/i / dy 2 R 1 {y ll y 2 ) 



I to IT Jto/T 



T ) \ ' 3 d / 2 + 1 
which can be exponentiated to give 



in 



1 f 1 

~ / dj/i / dy 2 R 1 (y 1 ,y 2 ) 

Z Jto/t! Jto/tl 



where 



-1/ *./ 

Z Jt /t2 Jt /t 2 



*(*1,*2) 



^2 



(! + •••) 



- CJ 



3 d/2+l 



(176) 



(177) 



dy 2 Ri{yi,y 2 ) 



1 rti/T ,i 2 /T 

*2 = ~ / cfe/i / dy 2 R 2 (yi,y 2 ) 

D Jto/T Jto/T 



- "4 / d 2/i / dy 2 R 2 {y x ,y 2 ) 

11 Jto/t! Jto/t! 

- To f dyx [ dy 2 R 2 {y x ,y 2 ) 

LZ Jt /t 2 Jt /t 2 



It is just this expression for A discussed in Section II. 

Turning next to f(x,t\/t 2 ), one of the first things to 
note is that it is for small x. Looking at the first term in 
the power series expansion in x 2 we obtain 



W{x,h/t 2 ) = -^x 2 



1 



l{ 2 \h/T) 



+ il[ 2 \t 2 /T) + hi 2 \ tl /t 2 ) 



and 



where 
/ 



I rti/T ,1/! 

W(x,t 1 /t 2 ) = - I dyx / dy 2 R x (yi,y 2 ) 



(2) f tl/T f Vl 

?>(ti/T)= dyx dy2Ri(yi,y2)Ag 1 (y 1 ,y 2 ) 

Jo Jo 



e 2 



Agix 2 _ ^ 



4 2 \h/t 2 ) = 



rti/T j-t 2 /T 

/ d Vl / 

Jo Jo 



dy 2 R 2 {yi,y 2 )Ag 2 (y 1 ,y 2 ) 



I rt 2 /T pyx 

+ 2 J dyi J dy 2 Ri{yi,y 2 ) 



I r^/T rt 2 /T 

g / d yi dy 2 R 2 (yi,y 2 ) 



e -±A gi x 2 _ 1 



\Ag 2 x z _ 1 



where since T ^> to, the lower limits can be set to zero. 
After making the change of variables y 2 — y\z one finds 
that one can perform the integral over y\ to obtain 



Ii 2 \t 1 /T) = -[J d (ti/T)-K d (t 1 /T)] 



where 



where 



A</i = 51(2/1,2/2) - l = 



3yi - y 2 - (yi - y 2 ) 2 



rti/T 

J d (h/T) = / dz 
Jo 



y d/2-l 



[(l + z) 2 (2-z)] d / 2 



and 



Ag 2 = g 2 (yi,y 2 ) - l = 



2/1 + 2/2 

3-2/1-2/2 



K d (h/T) 



ti/T 



dz 



yd/2-1 



Let us first look first at A<f>(ti,t 2 ) and the integrals 
$1 and $2. An investigation of $1 shows that it is a 



I< 2 > (ii/T) = 



[(l + z)(3-z)] d / 2 



£d(ti/t 2 ) + id(ii A2) - 3572 Afd 



(178) 



(179) 



(180) 



(181) 
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and 

L d {h/t 2 ) 



tl/*2 



dz 



r d/2-l 



where AVt^ 1 ) (x, t\/t 2 ) has no singularity for large x. If 
we now make the change of variables 



[(1 + ^)2(3-^(1 + ^/ 



(182) 



Notice that at equal times, Kd(l) = Kj, Ld(l) = Jd(l), 
and all of these integrals are well behaved. 

There is one last regularization that must be carried 
out before the perturbation theory expression for the cor- 
relation function can be used for all values of the parame- 
ters. Consider that W can be written as the sum of three 
terms: 



W^(x,t 1 /t 2 ) = ^J^ d m dy 2 R 1 (y 



,2/2) 



x e 



-±A gi x 2 _ 1 



W^(x,t 2 /t 1 ) = ^JJ / d yi J^ dy 2 Ri(yi,y2) 



x e 



r t 2 /T 



W<- 2 Hx 1 t 1 /t 2 ) = ^J^ / d yi J^ dy 2 R 2 {y u y 2 ) 



e 2 



±Ag 2 x 2 _ 1 



In W^(x, ti/t 2 ) let y 2 = y\z which leads to the result 

1 f U/T d yi 



W^{xM/h) = \J o 



yi 



f 

Jo 



x / dz 



yd/2-1 



[(l + z)(3-z- yi (l-z) 2 )]<V 2 



l x 2 «l(l-*)" 

2 X (3-^-y 1 (l-.)2) _ I 



(183) 



Notice in this integral that there is an apparent log di- 
vergence for small yi which is cancelled between the ex- 
ponential term and the subtraction term. Note that the 
log terms survives if x is large enough. To pick out the 
log term we can expand about small j/i except for the 
contribution x 2 yi and write 



ti/T 



wV{x,t 1 /t a ) = H- ^1 dz 



2 J yi Jo [(l + z)(3-z)]<V 2 



1 „2 

e -2 X (3-,) _ 1 



+ AW^{x,h/t 2 ) 



Vi 



Us 



T(l + x 2 ) 
in the leading integral we obtain 



(184) 



W^{xM/t2) = \J Q 



1 ; ds 
s 



dz 



yd/2-1 



[(l + z)(3-z)}^ 



1 x 2 «1 s(l-») 2 
g 2 i + x 2 T (3-z) 



1 



where the • • • refer to contributions which are regular. 
The integral over s can be divided into a regular part 
from to 1 and the singular part for large x given by 



yd/2-1 



z)(3-z)]<V 2 



1 x 2 tl a(l-z) 2 
g 2 T (3-z) _ I 



The singular part is now isolated in the second piece that 
does not have exponential convergence for large s . We 
have then 



W (1) (x,ii/*2) = ~ K d ln(l + x 2 ) + - 



Clearly 



^ (1) 0M 2 Ai) = -^ K d ln(l + x 2 ) + - 



(185) 



(186) 



Turning to the third contribution we again isolate the 
small yi and y 2 behaviors and expanding in x 2 in places 
which do not contribute to the singularity gives 



rt 2 /T 



dm 



d/2-1 d/2-1 1 

1 i/2 1 



(yi+V2) d 3 d ' 2 



e -^x 2 ( yi +y 2 ) _ ^ 



+ ■ ■ 



Next make the coordinate transformations 



yi 



2/2 



T(l + x 2 ) 



S2t 2 



(187) 



(188) 



T(l + x 2 ) 

which results in the leading contribution to the integral 

rl+x 2 r l+x 2 



ds 2 
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d/2-1 d/2-1 
S l S 2 

(s 1 + s 2 ) d 



x 2 ( s l f l I 5 2*2 \\ 



Again the leading behavior for large x comes from the 
terms which do not have exponential convergence for 
large s\ and s 2 : 



W {2 \xM/t2) = - 



1 



2 3 d/2+i J o 



l+x s 



l+x s 



ds-i 



d/2-1 d/2-1 
S-t So 

V —i ± L . . . 

(S1+S2Y 

If we let si — ► 1/si and s 2 —> l/s 2 then this integral 
takes the form 



W^(x,t 1 /t 2 ) = - 



1 



2 3 d / 2 +! J 1/{1+X l 



dsi 



ds? 



1/(1+^) 



X. THE ORDER PARAMETER CORRELATION 
FUNCTION 

A. Perturbation Expansion 

We turn next to the connection between the correlation 
function for the auxiliary field m and the order parameter 
correlation function. If we look at the problem using the 
transformation given by Eq.(^8|), we have 

C(12) = 0A(1)V(2)) 



= (*(1)*(2)) + (a(l)u(2)) + (u(l)<7(2)) + (u(l)u(2)) . 

The key point is that since u(l) vanishes exponentially for 
large |m(l)|, the averages over these fields are down by a 
factor of L~ 2 relative to the averages over the field ct(1). 
Thus (ct(1)o-(2)) « 0(1) as L(t) -► oo, while (<r(l)u(2)) 
and (u(l)cr(2)) are of (D(L~ 2 ) and (u(l)u(2)) of 0(L~ 4 ). 
In the scaling regime we have 



d/2-1 d/2-1 
Si So 
v —i ± L . . . 

(si + s 2 ) d 

This integral was evaluated in our treatment of with 
the result 



C(12) = (a(l)a(2)) . 



(192) 



W^(x, tl /t 2 ) = --^—2M d lr 



1 + x 



r-2 



2 3 d / 2 + 

Then to leading order for large x we have 

1 



W(x, h/t 2 ) = -K d ln(l + x 1 ) 



3 d/2+ 



T M d In 



The scaled correlation function then has the form 



f(x,t 1 /t 2 ) = f (x) 



l-u/2 



2 d+l 



M d 

3 d/2+l 



This can be exponentiated to obtain 
fo(x) 



f{x,h/t 2 ) 



[1 



(1 + a: 2 ) 
(189) 

(190) 



This quantity can evaluated using the two-point proba- 
bility distribution 

P h ( Xl x 2l 1, 2) =< S( Xl - m(l))5(x 2 - m(2)) > fc (193) 

and the correlation functions are obtained at zero exter- 
nal field via: 

C(12) = J dx x J dx 2 a{x 1 )a{x 2 )P {x 1 x 2l \2) . (194) 

More generally we can treat the set of correlation func- 
tions 

CW(12) = (ov,(l)0*(2)) 



dxi / dx 2 a n (xi)ai(x 2 )P (xix 2 , 12) 



where the exponent governing the large x behavior is 
given by 



Since we have computed the auxiliary field correlation 
functions to 0(2) we also need to determine C n g(12) to 
second order. This expansion can be developed as fol- 
lows. As in the case of the one-point quantity Ph(x±, 1), 
we again use the integral representation for the 5-function 
to obtain 



= c^ 2 2 d+1 



K d 



M d 



3 d / 2 



P h {x 1 ,x 2 ,U) = I ^ / ^ e -*-e-*- 



(191) 



2tt 



2tt 



< e 



H{12) 



>h 



where Ti(\2) = ifeim(l) + ik 2 m(2) . The average of 
the exponential is precisely of the form which can be 
rewritten in terms of cumulants: 



< e 



H{12 



> >h= exp 



£^4 n) (i2) 



(195) 
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where G&(12) is the n th order cumulant for the field 
H(12). These cumulants can all be expressed in terms 
of the m-field cumulants. In this section we can work 
directly in terms of zero external fields where all odd m- 
ficld cumulants vanish. We will need here 

G^ } (12) = (ifc 1 ) 2 G 2 (H) + (ik 2 ) 2 G 2 (22) + 2{ik 1 ){ik 2 )G 2 {12) 



d 2 d 2 

+ 6G 4 (H22) 2 2 + 4G 4 (1222) 

CLJL ^ (JLJL 2 



_d_d*_ 

dx\ dx\ 



P^( Xl ,x 2 ;12) 



(198) 



G^ } (12) = (i/ci) 4 G 4 (llll) + 4(^ 1 ) 3 (^ 1 )G 4 (1112) 
+ 6(i/ci) 2 (ifc 2 ) 2 G 4 (1122) + 4(i/ci)(ifc 2 ) 3 G 4 (1222) 
+ (ifc 2 ) 4 G 4 (2222) 

and 

G^ } (12) = (ifc 1 ) 6 G6(llllll) + 6(i/ci) 5 (^ 2 )G 6 (111112) 

+ 15(^i) 4 (i/c 2 ) 2 G 6 (111122) 
+ 20(^i) 3 (i/c 2 ) 3 G 6 (111222) + 15(ifci) 2 (ifc 2 ) 4 G 6 (112222) 

+ 6(ifci)(ifc 2 ) 5 G 6 (122222) + {ik 2 ) 6 G 6 (222222) . 

Working to second order in the (^-expansion requires 
keeping terms 



<e^ > = $ o( fc 1 ,fc 2 ,12)(l + -G^ + ^(-G 



1 



(4) 



1 ( 1 



4! 



2 I 4! 



+ 6! G « 



(196) 



where in zero external field $o{ki, k 2 , 12) is given by 
$o(fcifc 2 ,12)= e -^ fc ' G2 ( n ) +fe ' G2 ( 22 ) +2fclfc2G2 ( 12 )] . (197) 



where Pq°^(xi,2; 2 ; 12) is given by 



P^( Xl x 2 ,12) = 



1 



7(12) 



x e 



27r y/S 2 (l)S 2 (2) 

t2(12) , [^52(2)+x^52(l)-2G 2 (12)xix 2 ] 



"2S 2 (1)S 2 (2) I 



and 



7 (12) - 



/(12) = 



1 



V 1 - / 2 ( 12 ) 
G 2 (12) 



y/S 2 (l)S 2 (2) 



(199) 
(200) 

(201) 



Notice that we already have one resummation here since 
it is the full G 2 which appears in P^\x\x 2 , 12). The 
general set of two-point order-parameter correlation func- 
tions are given, up to 0(1) by 



G„,(12) = G^ ) (12) + Gi i ;(12). 



(202) 



where the first-order correction to the leading gaussian 
behavior is given explicitly by 



clV(12) = ffGr(llll)G^(12) 



1 



(o) 



+ -Gi 1) (1112)G^ £+1 (12) 
+ iGi 1) (1122)G^ 2! , +2 (12) + iGi 1) (1222)G^ M+3 (12) 



B. Terms of 0(1) 

Let us look first at the theory keeping terms up to 
O(l). It is then straighforward to show that the two- 
point probability distribution is given by 



^Gi 1) (2222)G^ +4 (12) 



(0) 



P {x u x 2 -12) = 



where we remember that G^(12) is a known functional 
of the exact G 2 (12) and the G 4 's are also the exact quan- 
tities. Thus G^ } (12) and G^?(12) both contain contri- 
butions of C(2). 

Thus determining the C n e to first order requires 
first determining G 2 to first order, then evaluating the 
G^(12) as functions of G 2 , and finally the evaluation of 
the various two-point contractions of G 4 (1234) evaluated 

at first order. We have already evaluated G 2 °^ and we 
have the appropriate first-order expression for G 4 (1234). 
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Notice that it is the same set of contracted four-point is the equal-time limit of Ldktiji'i) defined by Eq.(182) 
quantities which enter into the determination of any set Next we have 
n, I. Let us restrict our subsequent analysis to the cor- 
relation functions Coo (12). The first-order correction to 
the leading gaussian result 



Gf(1122) = -2w2 d S 2 (l)S 2 (2) 



C$(12) = -^«n-7(12) 



(203) 



(212) 



is given by 

c#(i2) = ^Gf(ini)c£to 



where 



..d/2-1 



[(t + y)(l+{ + y(l-y))] 



d/2 



i Gi 1) (1122)C^(12) + -G^^(1222)C^(12) 



Ml 



t-2 



-Gi 1) (2222)G^(12) 



(204) 



The zeroth-order correlation functions G^°](12) can be 
evaluated using the identities 



and 



dS 2 (l) 



^S(12) 



(205) 



(206) 



derived in TUG. Starting with the expression for Cq°q (12) 
given by Eq. ( [203|) , all of the other quantities can be cal- 
culated by taking derivatives. A summary of the results 



we need for Gq 1 q(12) is given by 

2Jl_ 
nSl(l) 



C|?(12) 



7 3 /(3-2/ 2 ) 



c<°>(i2) 



0o 2 



^S 2 3/2 (l)S 2 1/2 (2) 



7r5 2 (l)5 2 (2) 



7 3 / 



(207) 



(208) 



(209) 



Next we need the va rious contracted four-point cumu- 
lants appe arin g in Eq.(204). These can all be evaluated 
using Eq/( |15S| ). First we have the fully contracted four- 
point cumulant: 



G ( f\llll) = -4cj2 d S 2 (l)L d 



where 



Ld= dy 



y 



d/2-1 



o [(2-y)(l+y)2] 



2] d / 2 



(210) 



(211) 



x{ = r7(4fO 



and 



i + f + y(i-y) 

Finally we need 

Gi 0) (1112) = -w2%(l)5 2 (2)^ /4 - 1/2 

x [3Wi(xi,0 + W 2 (si.O] 



(213) 



d/2 



(1 + y)(l + 3£ + y(3-0 -22/ 2 ). 



x e -|*iff2(€,y) 



(214) 



Wa(xi,0= f dyy d ' 2 - 1 
Jo 



1 d/2 



(1 + y) 2 (l + 3£ - 2y) 



(i 



1 + 3C- 2tf) 



and 



4(1 + 



l + 3e + y(3-0-2y 2 ) 



Note t he c heck on Eqs.(212) and ( |213| ) that they reduce 
to Eq.( |21C| ) as 2 — > 1. Pulling all of these results together 
we have the 0(1) corrections to the order-parameter cor- 
relation function 



-2i d /(3 - 2/ 2 ) 
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+ £, d/A [3W 1 (x 1 ,0 + W 2 (x 1 ,0\ 

+ r d / 4 [3 w 1 (x 2 , r 1 ) + w 2 (x 2 , r 1 )] 



3c d/2 /w3(^i, o - 3r d/2 /^ 3 (x 2 , r 1 ; 



(215) 



There are several limits of interest. First consider the on- 
site unequal-time t\ ^> t 2 limit which gives the nonequi- 
librium exponent A. In this limit / is small and has the 
form for large £ given by Eq.(145) which, in our notation 
here, reads: 



(216) 



where — 2 d ^ 2 at lowest order. Similarly the zeroth- 
order order parameter correlation function has the form 



c(°)(o,t 1 ,i 2 ) = -^4 0) r d/4 



(217) 



Second-order contributions to / change the exponent A 
from d/2 to the expression given by Eq.(|l3|). Do the new 
first-order terms give contributions which change A at 
first order? To answer this question we need to set r = 
and take £ large in C^(12). The main results we need 
for large £ are: 

w 1 {o,t) = c d/2 <nW 



W 2 (0,0=C d/2 qo(d) 
W 1 (0,C 1 )=q2(d) 



W 2 {Q 1 C 1 ) = 1 J d d,2 C d ' 2 
a 



w 3 (o,o = -^r d 



W 3 (0,C 1 ) = [lnt + q 3 (d)} 
with the dimensionality dependent quantities defined by 



Qo(d) 



qi (d)= / dyy*' 2 - 1 
Jo 



q 2 (d) = / dy y^ 2 - 1 
Jo 



d/2 r2 



T{d) 



(i + y)(3-y) 



d/2 



2 d / 2 K, 



d/2 



(l + y)(l + 3y-2y 2 ) 



We see then that these contributions do not contribute 
terms to the exponent A at first order. One has only a 
contribution to the amplitude 



(1) (O,*i,i 2 ) 



lo2 c 



r d/4 [-6/oid + 3gi(d) + q (d) + 3q 2 (d)} 



For equal times t\ = t 2 = i, £ = 1, the expression for 
Cqq' ) (12) simplifies significantly. All of the W's share the 
same integrand except for the g's multiplying x 2 in the 
argument of the exponential. After considerable manip- 
ulation we have 



CS(12) = ^ 
on 



dy 



y 



d/2-1 



[(2-y)(l + y) 2 f 2 



(218) 



-2e"^ 2 (3 - 2/ 2 ) + 2e - 3* a M 1 .iO + 6(1 - f) e -h* 3 Mhv) 



For large x one has, since g (l,y) > 1, gi(l,y) > 1, that 
only the fully-contracted four-point cumulant contributes 
and 



C$(12) = ^(-6£ rf )/(12) 



(219) 



which does not give a first-order contribution to the ex- 
ponent v. 

The last point to be discussed in the 0(1) evaluation of 
the order-parameter correlation function is that one must 
be careful about the behavior at short-scaled distances at 
equal times. Let us write 



C-$(12) = 7 a A0r) 



(220) 



where A(x) can be read off from Eq.( |218| ). Then for 
small x we see that A (a;) goes to zero as x 2 while j 3 
blows up as 1/x 3 . Such singularities are unphysical and 
indicate that we are expanding about a singular order- 
parameter interfacial profile in this regime. It is just 
these singularities which give rise to Porod's law. One 
possible resummation is 



T{x) — —sin 

7T 



2 g 



— A(x) — 



f{x) 2 + 4q Q /ir 



(221) 



and a reasonable choice for go is 9o = ^(a:), where x is 
some renormalization point chosen such that the scaling 
function is smooth. 

We can then conclude that the first-order correction 
to the order parameter correlation function do not lead 
to corrections to the exponents A and v. Thus to 0(1) 
G 2 (12) and 0(12) share the same OJK exponents. 
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C. Terms of 0(2) 

Turning to the more involved case of the terms of 
0(2) contributions to C(12), we find four types of terms. 
There are the terms in G2(12) which are of 0(2) which 
must be included in the contribution given by Eq.(202) 
and which have already been evaluated. Inserting these 
results for / into Eq. ( 203 ) we find, since / is small in the 
regime controlled by v and A that the order-parameter 
correlation function picks up the same corrections leading 
to A m and v m . 

Let us turn now to the other three co ntrib utions: 

i). One has terms of 0(2) from Eq.(204) where the 



leading order G4 contributions are replaced by their 
next-order corrections G 



(2) 



ii). The leading Gq contribution in Eq.(196) is of 0(2) 
and contributes to 0( 2 )(12). 



iii). The term ( G^ 3 ) in Eq.(|l96|) is also of 0(2). 



Just as in Eq.(198) the factors of ki in Eq.(196) just 
lead to higher-order subscripts in the multiplicative fac- 
tors of C, ( l °](12). Thus in G (2) (12) there are contributions 
of the form 



dx\ <j{x\) 



f f dk 1 dk 2 _ 

J dX2 a[X2) J 
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x — (ifc 1 ) 3 (ifc 2 ) 3 G 6 (111222) 
6! 



order-parameter scaling function as in Ref.( ||). Let us 
return to the equation of motion for the order parameter 
given by Eq. (|39|) . We will limit the discussion for sim- 
plicity to the case of equal times where the equation of 
motion for the order-parameter correlation function can 
be written 

^ - 2V 2 R ) C(R,t) = -2< a 2 (l)(Vm(l))V(2) > 



-2JT(12) 



(222) 



Assuming we have a scaling solution, C(R, t) = F(x), we 
easily find that the equation of motion takes the form 



x ■ VF(x) + V 2 F{x) = L 2 K{\2) 



(223) 



We can then focus on if (12). At leading order it is easy 
to show that K(12) can be written in terms of the two- 
point probability distribution as 

K{12) =< (Vm(l)) 2 > G 20 (12) + K^(12) (224) 



whe re G 2 p(12) can be evaluated using Eq.(g03|) in 
Eq.(p|: 



C$(12) 



At lowest order we can also evaluate 



f(x)j(x) . 



(225) 



and 



G$(12)^G 6 (111222) 
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J dxt a( Xl ) J dx 2 a(x 2 ) J ^^ e -*-e-*- 



k {a \\2) = 



v (3) v (4) 



dxi / dx 2 <r 2 (xi)a(x 2 ) 



S 2 



5h{3)Sh(4) 



P h { Xl x 2 ,l2) 



|3=4=l,h=0 



- i I G^(1122) 



1 (Q{iki) 2 {ik 2 y 
2 



4! 



Ci°i(12) 



30 



2(4 



\\2 



G 2 (1122) 



While there are a great many terms to be analyzed from 
all three sets of contributions, none appear to give cor- 
rections to the exponents A and v. Thus it appears to 
0(2) that A = A m and v = v m . 



V (3) V (4) 



G 2 0) (13)G( 0) (14)G 4 o(12) 



+ G 2 0) (23)Gf (24)G 22 (12) 



+ (G 2 0) (13)G 2 0) (24) + Gf (23)G 2 0) (14))G 32 (12) 



|3=4=1- 



XI. EQUATION OF MOTION CONSIDERATIONS 

In the last section we determined the order-parameter 
correlation function by relating it to the auxiliary field 
correlation function. It is also instructive to consider 
using equation of motion methods for determining the 



Since 

(vfGf(13)) | 3 ^ = (v! 4) Gf(14)) | 4=1 =0 
we find 

X(°)(12) - (vf ) G 2 0) (12)) 2 G 22 (12) . 



(220) 
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Putting these results together, using Eq.([209| 



XII. CONCLUSIONS 



and 



f(xh(x) = tan[-F(x)] 



V x F(x) = -j(x)V x f(x) , 



(227) 



(228) 



we obtain the scaling equation 

f'VF(i)+V 2 F(3 



7T 

tan[—F(x) 



l-^(V x F(x)f 
[i 2 



= 



where 



L 2 < (Vm) 2 > 



2ft 



S. 



(0) 



(229) 



(230) 



Analyzing Eq. ( |229| ) , however, one finds the remarkable 
result that there is an analytical solution given by 



F = —sin 

7T 



(231) 



with \i given by the OJK result A* = 53 ■ Thus the 
equation of motion method gives the same result as the 
method of evaluating C^,(12) directly. This lends strong 
support to the structure of the approach. 

One interesting question concerns the work here com- 
pared to that in TUG. If one makes the replacement 



1 



1 



^-f(v^) 2 -- 

fl 2 fJL 



(232) 



in Eq.(229) one obtains the basic equation in the TUG 
approach where u must then be determined as the so- 
lution to a nonlinear eigenvalue problem. Given the se- 
lected u* one then has the analytic results: A = d — j 2 ^, 
and v = d — g^r . In the TUG approach A and v are 
not independent (y = 2\ — d). The TUG results for A 
and v as functions of d are shown in Tables I and II. It 
seems clear that the TUG values for A and v are superior 
to those of the present second-order theory. In particular 
the TUG approach gives the known exact results in one 
dimension. Clearly the behavior of v with d is different 
in the two approximations with TUG giving v — > and 
the current theory giving v — > 00 as d increases. Thus 
the TUG results are in agreement with the speculation of 
Bray and Humayun that large d corresponds to the OJK 
limit. 

It should be kept in mind that while the TUG ap- 
proach does well in giving these exponents it is deficient 
in treating the smoothness of the auxiliary field correla- 
tion function which enters into the determination of de- 
fect dynamics. The present theory gives good results for 
the exponents-the corrections to OJK all appear to be in 
the right direction and give the auxiliary field correlation 
functions which are smooth. 



The key accomplishment in this paper is to show how 
one can set up a systematic calculation which allows for 
both nontrivial noncquilibrium exponents, A and v 1 and 
results for the auxiliary field correlation function which 
are smooth enough to offer physical results for defect 
structures. The structure of the theory presented here is 
appealing since the deviation from the OJK results can 
be handled perturbatively and one can see and control 
the post-gaussian corrections. 

The key question remaining is the uniqueness of the 
theory developed here. While it seems that there is a 
degree of universality in this problem, the answer to the 
question posed involves the robust nature of this uni- 
versality. This can be investigated by looking at those 
changes in the equation of motion for the auxiliary field 
which may change scaling results. Thus one can try to 
find marginal variables which could lead to exponents 
which depend on a parameter. It seems reasonable that 
the uniqueness of the particular realization of the theory 
presented here can be investigated within its own struc- 
ture. 

It is clear that one can introduce more sophisticated 
resummation methods than the direct method used here. 
These methods could be useful, for example, in estab- 
lishing the connection between the exponents governing 
the statistics of the auxiliary field and those governing 
the statistics of the order-parameter field. While these 
are equal at 0(2), it is not at all obvious that this holds 
at higher order. It will be interesting to see if one can 
formulate alternatives to and resummations of the 4>p ex- 
pansion developed here. For example it seems to desir- 
able to find an expansion, like in TUG, which matches 
onto the exact solution for d = 1 . 

In the second paper in this series the theory is extended 
to the n- vector model. It will turn out, at least at 0(2), 
that the ^-expansion is related to a large-n expansion. 
The main focus of this second paper will be to look at the 
determination of defect spatial and velocity correlations. 
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APPENDIX: SOLUTION FOR S< 0) (T) 

In this appendix we present the solution to the equa- 
tion 

S 2 W-exp^^^^J Qd/2(t) 



where 



Q(t) = l + a(t- t ) 
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Cross multiply by Q d / 2 (t) and take the time derivative 
to obtain the rather simple equation 

6(0) (f ) | da Sj°\t) _ 

Introducing the integrating factor 

rt 



where 



(*) 



This is easily integrated to obtain 

X(t) = Si a) (t ) + 2Q jf* dt exp U j* dr-^ 

We can then do the integral 

/ dT 7^-r =01 ( dr— ln[l + a(r - t )] 
J t a Q( T ) J to dT 



= ln(l + a(t - to)) . 

This result is then inserted in the integrating factor and 
the remaining integral over i c an ea sily be carried out to 
obtain the result given by Eq.(134). 
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dimension 


theory 


OJK 


TUG 


best 


1 


0.6268.. 


0.5 


1 


l a 


2 


1.1051.. 


1 


1.2887 


1.246 ±0.02 b 


3 


1.5824.. 


1.5 


1.6726 


1.838 ±0.2 b 


large 


d/2 


d/2 


d/2 


d/2 c 



TABLE I. Values of exponent A from the current theory, 
OJK and TUG. a. Exact, see ref.( |l|), b. numerical results 
from ref.( c. Best guess. 



dimension 


theory 


OJK 


TUG 


LO 


1 


1.1596.. 





1 


0.4601.. 


2 


1.2492.. 





0.5774.. 


0.6877.. 


3 


1.3732.. 





0.3452.. 


0.9067.. 


large 


d(3 - 2^2) 








KV2-1) 



TABLE II. Values of exponent v from the current theory, 
OJK and TUG. The last column gives the values of the quan- 
tity UJ. 



